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Abstract 

Based on the numeric solution of a system of coupled channels for vector mesons 
(S- and D-waves mixing) and for tensor mesons (P- and F-waves mixing) mass 
spectrum and wave functions of a family of vector mesons qq in triplet states are 
obtained. The calculations are performed using a well known Cornell potential 
with a mixed Lorentz-structure of the confinement term. The spin-dependent part 
of the potential is taken from the Breit-Fermi approach. The effect of singular 
terms of potential is considered in the framework of the perturbation theory and 
by a configuration interaction approach (CIA), modified for a system of coupled 
equations. It is shown that even a small contribution of the -D-wave to be very 
important at the calculation of certain characteristics of the meson states. 



* e-mail : haysak@univ.uzhgorod.uzhgorod.ua 

* e-mail : feketeyu@univ.uzhgorod.ua * e-mail : morv@univ.uzhgorod.ua 

* e-mail : chalupka@pjsu.sk 



1 Introduction 



Meson states, which are considered as the bound states of a quark- ant iquark system, are 
convenient objects for studying both the strong interaction effects and various character- 
istics of weak interaction [1,2]. For the description of the low-energy properties of hadrons 
the following approaches are used: Bete-Solpiter approach, lattice calculation technique 
and potential models. Each of these methods has its own advantages and shortcomings. 
The potential models are the simplest from the point of view of mathematics, which is 
essential for practical calculations. In the framework of the present model the averaged 
mass spectrum [1], spin effects [3,4] and the decay widths of heavy quarkonia [5] can be 
well described. As concerning the light-quark systems, the situation is rather contradic- 
tory. Using the same parameters, some of the effects (spin effects, decay, averaged mass 
spectrum) can be described [6], but not all the effects together. The reason for this is the 
fact the light-quark systems being explicitly relativistic, and relativistic potential models 
should be applied for them. 

When the spin effects are considered in potential models (in Breit-Fermi approach), 
singular terms of the interquark interaction potential of the form ^ and S(r) arise. This is 
a serious problem at the calculations of the meson characteristics. As a rule, in such case 
perturbation theory is used [7]. But it has certain shortcomings. The main disadvantage 
is that the theory supposes small interaction, but in hadron physics it is the perturbation 
of the order 30-500 MeV, what is comparable with the distance between the unperturbed 
energy levels, hence, the condition of small perturbation not being fulfilled. The variational 
method also very popular in hadron physics [8]. These methods were in fact the first, 
which used to investigate the characteristics of mesons, considered as two-quark systems. 

In the triplet state configurations with uncertain orbital moment L — J + 1, L — J — 1 
occur. For example, mixing of S- and D-waves occurs in the state 1~ and P- and F- mixing 
for the state 2 + (we use the spectroscopic notations J p where the total moment and the 
system parity are indicated). In most papers where the triplet states are considered (See 
Ref. 3 and references therein) the authors neglect the channel coupling or introduce an 
additional parameter - the mixing angle [9-11]. 

In papers [12] the mixed state 1~ for gg-systems was considered in the frame of coupled 
channel and mass spectra with lepton decay width was investigated. It was observed a 
small Z)-wave mixing, but this small contribution causes essential influence on the value 
of hyperfine splitting and decay width. 

In the present paper the effect of the singular terms of the interquark interaction po- 
tential for the quarkonium state 1~, described by a system of coupled equations, is studied 
and the comparison of perturbation methods and configuration interaction, modified for 
the coupled equation system, is performed. We followed Ref. [12], when choose the 
Lorentz structure and parameters of gg-interaction with the numerical solution scheme 
of the coupled differential equation. The calculation version of the configuration inter- 
action approach based on the Schroedinger equation is taken from Ref. 3, where the 
hyperfine splitting is considered. In such approach the mixing angle for S- and D-waves 
is determined by the system dynamics and does not require any additional experimental 
data. 

2 Breit-Fermi approach 

It is widely accepted that the interaction between two quarks (or heavy quark-antiquark) 
consist of a short range part describing the one-gluon exchange and a infinitely rising 
long-range part responsible for confinement of the quarks 



K,(r) = -— + f3r. (1) 
r 

Wilson loop techniques suggest that the confining potential should be purely scalar, 
but relativistic potential calculations which have been published [13-15] show a need for 
(some) vector confinement. Maintaining Vy(r) + Vg{r) unchanged we do allow a fraction 
of vector confinement [6,13,14]. Namely 

W(r) = -^+A/r, V s (r)=(3 s r {{3 V + f3 s = f3) (2) 
r 

The confining potential transforms as the Lorentz scalar and vector potential trans- 
forms as the time component of a four-vector potential. As we can see, the choice of 
Lorentz structure of potential for quark- ant iquark interaction is important model for 
study of spin effects [7,14-16]. 

We consider the Breit-Fermi Hamiltonian for case m\ = m 2 = m 



with spin-orbit term 



H = V - + V + H LS + H ss + H T (3) 
m 



where S = Si + S 2 is the total spin of bound state and L is the relative angular 
momentum of its constituents, the spin-spin term 



and the tensor term 



H ss = -^AVV(r)S!S 2 (5) 



„ _ 1 fldVy d 2 V s \ 



where 



(Si r)(S 2 r) 

r 2 3 



, 12 = 12 ( ^ 11 - -SxS 2 ] (7) 



For bound state constituents of spin Si = S 2 = |, the scalar product of their spin, 
SiS 2 , is given by SiS 2 = — | for singlet states (S = 0), and SiS 2 = | for triplet states 
(S=l). 

Taking into account (2), then (4), (5) and (6) yield 

V LS= -^—(3^ + 3^-Sl) (LS) 

lm q mq y V T T J 



r ss = (— - 2na s S(r)) (S^) 

6m q mq \ r J 



The physical states of qq system are determined by total momentum J and its projec- 
tion Mj, parity P and total spin S (Table 1), here we use spectroscopic notation 2S+1 Lj 
Table 1. States of two fermion system 





Singlet states 


Triplet states 




(S- 


= 0) 


(S- 


= 1) 


J/p 


+ 




+ 









% 






1 


'Pi 




"Pi 


3 ^ + 3 D x 


2 




l D 2 


3 P 2 + 3 F 2 


3 D 2 



It is seen from Table 1 that there are singlet and triplet states with definite orbital 
moment and triplet states with mixed orbital components (L = J± 1). For pure states 
wave functions are of the form (F(r)/r) yf LS where yf LS is the spin-orbital part of the 
wave function. 

The radial part of wave function F(r) for singlet states satisfies equation 



d 2 F 



+ [k 2 







where l v r = m 1 Vn 



dr 2 r ,Zj 
V = V v (r) + V s (r) - (3/4)V ss and k 2 



(8) 



mE. 



The radial function of the pure triplet states Pq and states with L = J obeys the same 
equation that singlet states (8), but with potentials 3 V = V v + V S + jVss — 2V LS — 4Vr 
and 3 V = V V + V S + \V S s - V LS + 2V T respectively 

The wave function for ground triplet state of qq system with negative parity (P — —1) 
is a mixture of states 3 Si and 3 D X and may be put in the form 



Then the equation (H — E)ip = is equivalent to coupled system [17] 



(9) 



1 d 2 

[ — - E + 3 V c \u + VSV T W = 0, 

m dr z 



where 



1 d 2 6 

[ — -E+ + 3 V c - 2V T - 3V LS }w + V8V T u = 0, 

m dr z mr z 



1 



(10) 



3 V C = V v + V s + -V ss . 



The Schroedinger equations are linked due to the presence of a tensor component 
Yr in the interaction potential. In Ref. 12 the system (10) was solved numerically for 
regular part of the potential (4-6), and irregular terms were calculated as a first order of 
the perturbation theory. Hereinafter we study the contribution of the singular terms in 
the framework of the perturbation theory and CIA. We generalize these methods to the 
system of equations. For this purpose that is convenient to rewrite the system (10) in a 
matrix form. The Hamiltonian of the system can be separated to regular and irregular 
parts H = H + W, namely 



j_ d 2 

m dr 2 



-^ + kr 



3m' 2 r ' 



m dr 2 



f + kr + 



+ 



3m 2 r mr 



3 iWv-f3s' 



(11) 



and 



W - ( 3m"^)' ^mVi | (-to) 

VV — I 7g 3a . /g 3a _ _6o 9q . I 

\ * m 2 r 3 ' " m 2 r 3 m 2 r 3 2m 2 r 3 ' / 

Then unperturbated eigenfunction of equation 

H ^ n = E n m n (13) 



where 

*• = ( i: ) (i4) 

can be used as a basis for the full wave function of the system, namely, 

^ = E«A (15) 
Then the matrix elements of the perturbated term 

W mn = J ** m W* n d t (16) 



are given by 

W mn = J u m (-^^8(~t))u n dr + J w m (y/8^^)u n dr + J u m (y/8^^)w n dr + (17) 

7 ^ 3m 2 ^ ' m 2 r 3 2m 2 r 3 ^ ~~ 55 5T 5T ^ss+ls+st 

and in the framework of the first order of perturbation theory the correction to the 
energy spectrum is equal to 

A£ nn = W^" " + 2W^T n " + W^5+ls + st (18) 

It can be noted that mean value of unperturbed Hamiltonian Hq (11) consist of the 
same parts determined by S-wave, S-D-interferention and D-wave correspondent 

= (i>\H \4>) = E S + E SD + E D . (18a) 

As seen from Eq. (17), the S'-wave gives pure contribution only for the spin-spin interac- 
tion (the first term), the interferential SD term (the second and third ones) contains the 
spin-tensor correction, and pure D-w&ve (the fourth term) includes all the spin-dependent 
components of the interaction potential. This structure is reflected in Eq. (18) by sub- 
and superscripts. In the case of the configuration interaction method we used an algo- 
rithm, developed in Ref. 3,4. The mean value E of Hamiltonian H + W is an eigenvalue 
of the following system of linear algebraic equations: 

ai (E - El - W u ) - a 2 W 12 - a 3 W 13 - ... - a n W ln = 
-aiW 2 i + a 2 (E - E° 2 - W 22 ) - a 3 W 23 - ... - a n W 2n = 

-aiWm - a 2 W n2 - a 3 W n3 - ... + a n (E - E° n - W nn ) = 

where E®- are eigenvalues of Hamiltonian H (11) and Wij - are correspondent ma- 
trix elements (16). Respectively, eigenvectors (ai, a 2 , a n ) gives eigenfunctions <3>j, % = 
1, 2, n of the Hamiltonian H + W. 



3 Hyperfine and fine splitting 



To study the role of tensor forces we calculated for triplet states L = J ± 1 in mixed 
form (10), and for single S-wave (as most authors do). It was pointed [14,16] that the best 
agreement with experimental data for potential (1) is obtained when f3 = 0.18GeV 2 . For 
description of gg-system we were varying parameter (3y (and fixing Pv + Ps = 0.18GeV 2 ) 
) to achieve the agreement with experimental mass splitting of 1 states. Finally 
we have used the following parameters: f3y = O.OOlGeV 2 , fts = 0.179GeV 2 ior_uu- 
systems; f3y = 0.04GeV 2 , fis = 0.14GeV 2 for charmonium and bottonium, a s (bb) = 
0.24, a s (cc) = 0.38, a s (uu) = 0.54. The quark masses are: m b = A.7GeV m c = lAGeV, 
and m u = 0.33GeV. Tables 2-4 list numerical results, namely, mass spectra for single 
S'-wave and mixed state, fraction of tensor part of the potential to energy levels (18a) 
and D-wave fraction in the total wave function (P D = J \ w\ 2 dr). And for comparison we 
showed experimental values [19] and calculation with screened potential [3,4]. 



Table 2. Hyperfine splitting for the charmonium 



State 


S'-wave 
E theor ,MeV 


S-D-waves 
E theor , MeV 


[3] 

E theor , MeV 


[19] 

E exp , MeV 


% 


Pd, 

% 


1% 


2980 






2980 






l 3 ^ 


3153 


3097 




3097 


16 


0.05 


l 3 ^ - 1% 


173 


117 


110 


117 






2% 


3642 






3590 








3759 


3734 




3685 


3 


0.8 


2 3 5*i — 2^5*0 


117 


92 


67 


95 






3 1 Sq 


4107 












3 3 S 1 


4208 


4192 




4040 


1 


1.3 


3^Si — 3^5*0 


101 


85 










Table 3. Hyperfine splitting for the bottonium 


State 


S-wave 
E theor , MeV 


S-D-waves 
E theor) MeV 


[3] 

E theor , MeV 


[19] 

E exp , MeV 


E SL> , 

% 


Pd, 

% 


I 1 S 


9415 












l 3 Si 


9462 


9460 




9460 


2.1 


0.004 


l 3 £i - I'So 


47 


45 


46 








2'So 


9883 












2'S 1 


9911 


9911 




10023 


0.2 


0.04 


2 3 5*i — 2^5*0 


28 


28 


26 








3^0 


10201 












S 3 ^ 


10224 


10224 




10355 


0.1 


0.1 


3 3 5*x — 3^5*0 


23 


23 











Table 4. Hyperfine splitting for (ml)-systems 



State 


,S-wave 
E theor: MeV 


SD-waves 
E theor , MeV 


[3] 

E theor ., MeV 


[19] 

E exp , MeV 


% 


Pd, 

% 




140 






140 






i s s l 


674 


640 




770 


4 


0.02 


l 3 ^ - l'So 


534 


500 


923 


630 






2^0 


1134 






1300 






2 3 ^ 


1564 


1543 




1450 


1 


0.04 


2 3 5i - 2 1 S , 


430 


409 


411 


150 







It can be noted that we obtain good description of mass spectra and hyperfine 
splitting for heavy quark systems. For light system our quasirelativistic model is less 
success, here relativistic kinematics play sufficient role [6]. 

In Tables 5-7 fine splitting data are listed for P-states, where we indicate in brackets 
result with mixed tensor states ( 3 P2 — 3 F2 mixing). Comparison results for S — D and 
P — F system show that influence of tensor forces decreasing when total moment of the 
system increases. 



Table 5. Fine splitting for cc— systems (in MeV) 



State 


AM 


Our results 


[7] 


[20] 


[21] 


[22] 


[19] 


IP 


M( 3 P 2 ) - M( 3 Pi) 


51(49) 


576 


51 


45 


28 


45.67 




M( 3 P X ) - M( 3 P ) 


72 


76 


83 


91 


32 


95.51 




M( 3 P 2 ) - M( 3 P ) 


123(121) 


132 


134 


137 


66 


141.18 



Table 6. Fine splitting for bb— systems (in MeV) 



State 


AM 


Our results 


[7] 


[20] 


[21] 


[21] 


[19] 


IP 


M( 3 P 2 ) - M( 3 P0 
M( 3 Pi) - M( 3 P ) 
M( 3 P 2 ) - M( 3 P ) 


12(11) 

15 
27(26) 


23 
26 
49 


31 
41 
72 


24 
37 
61 


28 
32 
66 


19.9 ± 1.1 
32.8 ± 1.5 
52.7 ± 1.5 


2P 


M( 3 P 2 ) - M( 3 Pi) 
M( 3 Pi) - M( 3 P ) 
M( 3 P 2 ) - M( 3 P ) 


10 (10) 
15 

25(25) 


16 
20 
36 


24 
32 
56 


17 
26 
44 


20 
24 
44 


13.3 ±0.9 
23.1 ± 1.1 

36.4 ± 1.0 


Table 7. Hyperfine splitting in P— waves (in MeV) 




State 


AM 


Our results 


[7] 


[19] 




cc(lP) 
bb(lP) 
bb(2P) 


M( 3 Pi) - M( x Pi) 
M( 3 Px) - M( l Pi) 
M( 3 Pi) - M( l Pi) 


-10 
-3.4 
-3.5 


13 
4.3 
3 


-15.63 ±0.36 



4 Vector mesons in CI approach 

We also analyzed influence of perturbation part of the system Hamiltonian because most 
authors study spin-spin interaction in the frame of perturbation theory of the first order. 
To estimate high order theory we use the configuration interaction approach (CIA), the 
one configuration case of which coincides with the first order perturbation theory. 

Below the numerical results for the mass spectrum of cc and uu systems are given. The 
results of the mass spectrum calculations listed in Table 8 (nonperturbated eigenvalues 
of energy, correction according to the perturbation theory, energies according to CIA 
with two and three configurations). It is seen from the table that the series of Eq. (10) 
converges rapidly for heavy systems and already the perturbation theory gives more than 
90% contribution of the singular terms, but for light mesons the extension beyond the 
perturbation theory is essential. 



Table 8. Vector meson mass spectrum 



States 


M TH , MeV 


M TH ,MeV 


M TH ,MeV 


M TH ,MeV 


M EXP ,MeV 




non-perturb. 


1 conflg. 


2 conflg. 


3 conflg. 




J/M1S) 


3043.45 


3098.58 


3097.11 


3096.87 


3096.87 


i>(2S) 


3649.04 


3674.02 


3675.50 


3674.87 


3685.96 


^(4040) 


4098.04 


4113.11 




4113.97 


4040 


p(770) 


685.5 


778 


771 


768.5 


768.5 


p(1450) 


1574.5 


1643 


1650 


1643 


1465 


p(1700) 


2270.5 


2332 




2341 


1700 



Table 9. Perturbation theory for a J/ip meson 





W£§ Un ,MeV 


2W^ u ",MeV 


WSJSr +LS ,MeV 


AEn,MeV 


6.04 


49.6 


0.00050 


AE 22 , MeV 


4.18 


21.16 


0.00036 


AE 33 , MeV 


3.65 


11.95 


0.00052 



Table 10. Perturbation theory for a p meson 





W£g Un ,MeV 


2W^ u ",MeV 


WSSs T+LS ,MfiV 


AEn, MeV 


60.98 


31.67 


-0.19 


AE 22 , MeV 


50.19 


19.02 


-0.28 


AE33, MeV 


46.64 


14.69 


-0.075 



The absolute values of corrections due to certain singular terms of the potential in the 
framework of the perturbation theory are rather interesting. The correction values for the 
spin-spin, spin-tensor and spin-orbit components of the singular part of the interaction 
for charmonium and p meson are given in Tables 9 and 10, respectively. It is seen that for 
heavy systems the spin-tensor part of the correction value by order of magnitude exceeds 
the spin-spin part and is totally determined by the presence of the D wave admixture. For 
light systems the spin-spin and spin-tensor correction values are of the same order. Hence, 
when the spin effects are considered, one should take into account the orbital structure 
of the meson states, especially for the systems of heavy quarks. However, we note that 
the .D-wave admixture for the light and heavy systems is less than 1 % and about 4% , 
respectively [12]. In comparison with Ref. 10 , we note that for the ty(2S + D) state the 
D-wave admixture in our case is Pd = 0.008 [16], which corresponds to the mixing angle 
of (p = 5°. In Ref. 7 the mixing angle of pure triplet states 2S and ID (p = 12° is quoted. 



5 Leptonic decay of heavy quarkonia and wave func- 
tion in the origin 

For the leptonic decay widths of two-quark system we shall consider decay of vector 
states into e + e~ pairs. The leptonic decay width of system M qq - — > e + e~ is calculated 
from the Van Royen-Weisskopf formula [23] 

Ao 2 O 2 

r( 3 Si-e + O = -^|#(0)| 2 , (20) 

11 

where M qq - is mass of vector meson, Q is quark charge, a em is the fine structure con- 
stant and R(0) is the radial wave function in the origin. In our case radial wave function 
in the origin is determined by S'-wave component of the system wave function, namely, 
by ^ (9) because the D-wave component vanishes in the origin. 



The formula (20) is based on the notion that constituent quark-antiquark pair anni- 
hilates into a single virtual photon, which subsequently gives rise to a leptonic pair. The 
relativistic and radiative QCD corrections [24] modifies eq.(20) 



r( 3 Si 



e e 



f(l- 



16a, 



m„ 



3tt 



(21) 



As Eichten and Quigg have pointed out [25] the QCD correction reduces the 
magnitude of V significantly, however the value of reduction is somewhat uncertain . For 
vector mesons containing light quarks this formula leads to paradoxes [26]. In paper [27] 
Motyka and Zalewski extrapolated eq. (21) by rational and exponential function and 
obtain averaged formula. 



r 



V- 



[22) 



with F(c) = 4.73. 10 -5 for charmonium and F(b) = 2.33. 10~ 5 for bottonium. We have 
calculated decay widths using the formula of Van Royen-Weisskopf (20) and formula (22). 
Table 11 lists this results. 

Table 11. The leptonic decay widths of heavy mesons 



State 


S-waves 


SD-w&ve 


[27] 


[5] 


[27] 


[19] 




r theor .,keV 


r theor ,,keV 


r i/ie0 r.,keV 


T theor ,,keV 


r theor .,keV 


Texp. ;keV 


J/iplS 


8.2 (5.63) 


7.8 (5.41) 


4.5 


4.24 


8.0 


5.26 ±0.37 


i)'2S 


4.0 (2.79) 


3.7 (2.59) 


1.9 


1.81 


3.7 


2.12 ±0.18 


^"3S 


2.9 (2.01) 


2.6(1.82) 




1.22 




0.75 ±0.15 


TIS 


1.2 (1.01) 


1.14(0.96) 


1.36 


0.85 


1.7 


1.32 ±0.04 


T'2S 


0.63 (0.53) 


0.58(0.49) 


0.59 


0.38 


0.8 


0.52 ±0.03 


r"3S 


0.49 (0.42) 


0.44(0.37) 


0.40 


0.27 


0.6 


0.48 ±0.08 



The calculations of widths were performed with tensor forces and without ones. Value of 
the widths which were calculated by formula f!22l) are given in parentheses. 
Table 13. Wave functions in the origin |_R(0)| 2 . (In GeV 3 ) 



State 


\R(0W, 


\R(0)\\ 


\R(0W, 


\R{0)\*, 


\R{0)\», 


\R(0)\\ 


\R(0)[\ 




non-perturb. 


1 config. 


2 config. 


3 config. 


[5] 


[25] 


[18] 


j/^is) 


0.77 


0.69 


0.83 


.85 


1.05 


1.45 


0.81 




0.53 


0.56 


0.47 


0.50 


0.63 


0.93 


0.53 


^(4040) 


0.47 


0.53 




0.41 


0.52 


0.79 


0.46 


p(770) 


0.11 


0.086 


0.126 


0.136 








p(1450) 


0.091 


0.081 


0.073 


0.089 








p(1700) 


0.084 


0.12 




0.061 









The values of the squared radial wave function in the origin for the non-perturbed and the 
perturbed case along with the CIA calculations with two and three configuration sets are 
listed in Table 13. It is seen that, contrary to the energy spectra, the consideration of the 
decay widths beyond the framework of the perturbation theory results in the correction 
value of 20-30% . Our results are very close to those of Ref. 17 where the same quark 
systems were considered on the base of Schroedinger equation with the generalized Breit- 
Fermi potential. The difference from our work is in the choice of the Lorentz structure 
of the quark-quark interaction and neglecting the D-wave admixture. Also squared mean 
value radius obtained from our wave function coincide well with other calculations (Table 
12). 

Table 12. Squared mean value radius 



nL 


[22] 
(r 2 ) 1 / 2 , fm 


Our results 


[29] 
(r 2 ) 1 / 2 , fm 


Our results 
(r 2 ) 1 / 2 , fm 


1 SD 


0.43 


0.433 


0.24 


0.256 


2 SD 


0.85 


0.847 


0.51 


0.552 


3 SD 


1.18 


1.182 


0.73 


0.768 



It must be noted that squared mean value radius 0.7 fm meet the requirements of 
quark- ant iquark pair creation (uu) (string break). It is necessary to modify the potential 
model taking into account the opening of a new channel for those conditions where the 
given value is overlapped. 



6 Conclusions 

In the present work we have studied the influence of the quark-antiquark interaction 
potential structure to the meson properties. We pay the main attention to the role of the 
tensor forces. It was considered mass spectra and leptonic decay width. 

The regular part of the potential was taken into account by numerical solutions of 
Schroedinger equations for singlet and unmixed triplet states but a coupled differential 
equation for mixed triplet states (S — D- for vector mesons and P — F- for tensor mesons). 
The irregular part of potential was taken into account in the frame of the configuration 
interaction approach. 

The analysis of our results presented in Table 2-4 shows that differences between 
the theoretical calculations of mass spectrum for heavy quarkonia and the experimental 
results are 1-4%. Therefore we expect relativistic correction in the range 1 - 16 % for 
charmonium, and up to 1 % for bottonium spectrum. We have also been able to describe 
hyperfine splitting of cc- and 66-quarkonium. For describing mass spectrum of light mesons 
it is necessary to use relativistic potential model. We have calculated hyperfine splitting 
uu system, too. 

It shown that, admixture of D-wave is small (less then 1%) but its influence to the mass 
spectra reach 16%. The leptonic decay widths suggest that for charmonium theoretical 
widths, which were calculated by Van Royen-Weisskopf formula are systematically higher 
than experimental data. But for bottonium we have obtained lower values. For J / -^-meson 
better decay widths are obtained by means of formula ()21j) . which takes into account QCD 
correction. Here the influence of .D-waves ranges from 4 % for ground state up to 25 % 
for the second excited state, but for values calculated by Van Royen-Weisskopf formula it 
is from 8 % up to 50 %. For T-meson it is opposite. More exact is the calculation done 
by (j2Tj) . but QCD correction reduced the results more significantly. Besides, .D-waves 
contribute less than for charmonium: from 4 % for ground state up to 11 % for second 
excited state. 

The results obtained show that contribution of the .D-waves is impossible to neglect 
for considered leptonic decay width of quarkonia. 

It is shown that, in spite of a small admixture of the D-wave this component of the 
wave function plays an essential role at the account of the irregular part of the interaction 
potential, and namely the irregular terms are considered by most authors while the spin 
effects being studied. The presence of the D -wave essentially enhances the contribution 
of the spin-tensor component into the mass spectrum. As concerning the technique of the 
account of the irregular part of the interaction it should be noted that in case the mass 
spectrum of the systems being considered, one can restrict themselves by the first order 
correction of the perturbation theory, but at the analysis of the decay widths the CIA 
results essentially improve the first order perturbation theory. 
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